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Abstract. A precise measurement of the Dalitz plot parameter, a, for the rj — > 37r° decay is presented. 
The experiment was performed with the Crystal Ball and TAPS large acceptance photon detectors at 
the tagged photon beam facility of the MAMI-B electron accelerator in Mainz. High statistics of 1.8 ■ 10 6 
rj — > 37T° events were obtained, giving the result a = —0.032 ± 0.002 sta t ± 0.002 syst . 

PACS. 25.20.Lj Photoproduction reactions - 13.60.Le Meson production - 14.40.Aq n, K and rj mesons 
- 13.25.Jx Decays of other mesons 



1 Introduction 

The 77 — > 3ir° decay violates isospin symmetry. There- 
fore, it offers a unique possibility to study symmetries 
and symmetry-breaking characteristics of strong interac- 
tions. Because electromagnetic contributions to the am- 
plitude can be neglected jUGJE], this decay occurs due to 
the isospin breaking part of the QCD Hamiltonian: 



^( mu ~ ma) (mi - dd). 



(1) 
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Therefore, the amplitude is proportional to the mass dif- 
ference m u — md of the two lightest quarks u and d. 

Calculations of the decay amplitude are usually based 
on the framework of Chiral Perturbation Theory (%PT). 
The leading order term 0(p 2 ) of the momentum expan- 
sion yields the constant amplitude [4]: 

A(^ 37 r°) = g ° ( ^; 2 md) ~(m u -m d ), (2) 

where Bq and the pion decay constant F„ are the low- 
energy constants of xPT. There is further theoretical work 
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Table 1. Theoretical results for the Dalitz plot parameter, a. 



Calculation Refs. a 



XPT 0(p 2 ) 


KB 





XPT 0(p 4 ) 
XPT G(p 6 ) 


HI] 


0.015 




0.013 ± 0.032 


Dispersion 


7 


-0.007... -0.014 


UCHPT 


IS] 


-0.031 ± 0.003 



determining the amplitude in the second i5j and the third 
[6] chiral order, both of which include loop diagrams de- 
scribing hnal-state rescatterings of the pions. Higher or- 
der rescattering effects are examined by using dispersion 
methods [TUB] or the unitarised xPT approach (UCHPT) 
[9] based on the Bethe-Salpeter equation. 

The squared absolute value of the decay amplitude 
may be expanded around the centre of the Dalitz plot 
[TO] and Bose symmetry dictates the form: 

\A( V ^3n Q )\ 2 = \N\ 2 {1 + 2az + ...). (3) 

N is a normalisation constant, which is equal to the am- 
plitude that would apply, if the decay proceeded only ac- 
cording to the available phase space. The Dalitz plot pa- 
rameter, a, describes the pion energy dependence of the 
squared absolute value of the decay amplitude up to first 
order of the expansion. The parameter z is given by [lOj 

, =6 £(|Lz^y = _£_. ( 4) 

^\mr,-3m^oJ p l max 

Here Ei represents the pion energies in the rj rest frame 
and p is the distance from the centre to a point in the 
Dalitz plot. p m ax is the maximum value of p. z varies 
from z = 0, where all three pions have the same energy 
Ei = m v /<S, to z = 1, where one of the pions is at rest. So, 
determining the Dalitz plot parameter offers a nice pos- 
sibility to test the different theoretical calculations based 
on xPT (see tabled]). 

In lowest order of xPT there is no final-state inter- 
action between the three pions. Thus, the amplitude is 
constant [J] and a = [TJJ . With the amplitude including 
one-loop diagrams [5; Bijnens [11] calculated a — 0.015. 
Including two-loop contributions, Bijnens and Ghorbani 
[5] get a = 0.013±0.032. Contrary to the experimental re- 
sults summarised in tabled] all calculations based on xPT 
alone predict a positive value. Using instead dispersion re- 
lations [7] in combination with extended Khuri-Treiman 
equations [T2] gives the result —0.014 < a < —0.007, de- 
pending on the input parameters of the calculation. The 
unitarised xPT approach [9] based on the Bethe-Salpeter 
equation constrains its free parameters to known branch- 
ing ratios and ttt: scattering phases from [131 . The fit yields 
a = —0.031 ± 0.003, within one standard deviation in 
agreement with the high-statistics experimental results of 
the Crystal Ball at BNL [14 and KLOE collaborations 



Table 2. Experimental results for the Dalitz plot parameter, 
a. The result of Baglin et al. [16] has been omitted due to a 
lack of statistical significance. 



Experiment Refs. a 



Crystal Ball at BNL 


M 


-0.031 ± 0.004 


KLOE 


□S 


-0.027 ± 0.004lo;Soe 


GAMS-2000 


El 


-0.022 ± 0.023 


Crystal Barrel 


18 


-0.052 ±0.017 ±0.010 


SND 


19 


-0.010 ±0.021 ±0.010 


CELSIUS/WASA 


20 


-0.026 ±0.010 ±0.010 


WASA at COSY 


m 


-0.027 ± 0.008 ± 0.005 



Table [2] summarises the experimental world data set 
for the Dalitz plot parameter, a. The Crystal Ball collab- 
oration [14] measured a = —0.031 ± 0.004 from a sam- 
ple of 9.5 • 10 5 events obtained at BNL, while KLOE [15] 
found a = -0.027 ± 0.004+°;°^ with 6.5 • 10 5 events. All 
other experiments listed in table[2]collected not more than 
1.2 • 10 5 r] -> 3tt° events. The PDG [10] currently quotes 
a = —0.031 ± 0.004, but this value is dominated by the 
result of the Crystal Ball collaboration, since the latest 
publication of the KLOE collaboration has not been in- 
cluded. 

When proposing the measurement of the r\ — > 37r° de- 
cay at MAMI, only the result of one high-statistics ex- 
periment, the Crystal Ball at BNL, was published. This 
lack of precise data led to the conclusion that a new high- 
statistics measurement of the Dalitz plot parameter was 
needed, again using the Crystal Ball, but with the r\- 
photoproduction mechanism instead of a hadronic pion 
beam at BNL. This need for another precise measurement 
of the Dalitz plot parameter was emphasised even more, 
when the KLOE collaboration announced their prelimi- 
nary result a = -0.014 ± 0.004 stat ± 0.005 syst [22], which 
showed a big discrepancy with the BNL value. The KLOE 
collaboration at DA</>NE used e + c~ — > <p — > VI as the pro- 
duction reaction. Although the revised KLOE result [T5] 
agrees with the Crystal Ball at BNL value within the given 
errors, the fundamental importance of the Dalitz plot pa- 
rameter still requires further experimental input. 

This paper describes an analysis of tagged photon ex- 
periments carried out at the electron accelerator MAMI-B 
(Mainz Microtron B) [23U24j in the years 2004 and 2005 
within the Crystal Ball/TAPS collaboration [25] . Data 
from two different experiments with the CB/TAPS-setup, 
that differed in the trigger conditions and the tagged pho- 
ton energy range, were analysed. The first was dedicated 
to the neutral decays of the r\ meson and especially to the 
rare r] decays. The second investigated radiative tt° -photo- 
production. 



2 Experimental setup 

The Dalitz plot parameter, a, was determined from 37r° 
decays of r\ mesons produced in the 7p — > ?yp reaction. 
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Fig. 1. The Crystal Ball/TAPS setup. 



The photons were emitted by bremsstrahlung of 883 MeV 
electrons from the MAMI-B [23p4] accelerator. The elec- 
trons were separated from the photons and momentum 
analysed by the Glasgow tagged photon spectrometer [26l 
12 T] at Mainz, with an energy resolution of approximately 
2 MeV. The photon energies were determined by energy 
conservation with a maximum photon flux above the 77 
threshold of roughly 1 • 10 5 7/(s MeV). 

The 4.8 cm long liquid hydrogen target was located 
at the centre of the spherical Crystal Ball (CB) [281129] 
photon spectrometer. The CB consisted of 672 optically 
insulated Nal(Tl) crystals, each read out by an individ- 
ual photomultiplier. It covered the full azimuthal angle 
range for polar angles between 20° and 160°. Each crys- 
tal had the shape of a 41cm (15.7 radiation lengths) long 
truncated pyramid pointing towards the centre of the CB. 
Electromagnetic showers were measured with an energy 
resolution of a/E 7 = 0.02/(E 1 /GeV) 1 / 4 . The resolution 
in the polar angle ag was 2° to 3°, while for the azimuthal 
angle it was a^ — ag/ s'mO. The target at the centre of 
the CB was surrounded by a thin particle identification 
detector (PID) [30] to register charged particles hitting 
the CB. It consisted of 24 30 cm long and 2 mm thick op- 
tically isolated plastic scintillator bars arranged parallel to 
the beam axis, so that each covered 15° of the azimuthal 
angle range. 

The forward angles between = 4° and 8 = 20° were 
covered by TAPS [3TU32] consisting of 510 BaF2 crystals 
arranged as a wall, preceded by a single layer of 510 5 mm 
thick veto plastic scintillators. The CB and TAPS together 
covered roughly 97 % of the total solid angle (fig.Q]). TAPS 
was positioned 173 cm downstream of the centre of the 
CB, giving the opportunity for an efficient time-of-flight 
analysis for particle identification. Each of the hexagonally 
shaped BaF2 crystals had an inner diameter of 5.9 cm and 
a length of 25 cm, which corresponds to approximately 
12 radiation lengths. The shower energy resolution was 
a/E 1 = 0.0079/(S 7 /GeV) 5 + 0.018. The angular reso- 
lution for 300 MeV photons was 0.7° full width at half 
maximum (FWHM). 

The experiment trigger comprised two event condi- 
tions. First, the total sum of the CB photomultiplier ana- 
logue signals had to exceed a threshold that corresponded 



to approximately 390 MeV for one experiment and 60 MeV 
for the other. Secondly, the sector multiplicity in the CB 
and TAPS had to be greater than 2. Here up to 16 adja- 
cent crystals made up a sector in the CB, giving 45 hard- 
wired sectors. In TAPS each sector was one quarter of the 
wall. If at least one of the signals in the sector exceeded 
a threshold of 20 to 40 MeV, depending on the relative 
calibration of the photomultiplier signals, it contributed 
to the multiplicity. 

The determination of the Dalitz plot parameter re- 
quires the elimination of the phase space contribution to 
the amplitude. Therefore, a Monte Carlo simulation of the 
experiment was produced with a = 0. For the analysis 
described in this paper 100 • 10 6 rj — > 37r° events were gen- 
erated and the particle propagation simulated with the 
GEANT v3.21 software. These events were analysed in 
the same way as the experimental data, resulting in an 
average reconstruction efficiency of £,,^3^0 w 23 %. 



3 Data analysis 

The Dalitz plot parameter, a, for the decay 77 — > 37r° was 
calculated from the analysis of the reaction 

7P — * ?7p — * 37r°p — * 67P. (5) 

The first step of the event selection demanded six clus- 
ters in the CB and TAPS identified as photons, ignoring 
the other particle types like protons, charged pions and 
neutrons. With the CB and the PID, clusters were iden- 
tified as charged particles by checking the agreement of 
the azimuthal angles of CB clusters with the 4> angle of 
the hit PID elements. Photons had no corresponding PID 
hit. The charged particles were then divided into protons 
and charged pions by comparing the energy of the clusters 
with the energy deposit in the PID scintillator. Neutrons 
could not be separated from photons in the CB. In TAPS 
photons, protons, charged pions and neutrons could be 
identified using information from the veto detectors, com- 
paring the time of flight with the energy deposition in 
the BaF2 crystals and analysing the pulse-shape of the 
BaF 2 signals. A cluster was formed by a group of adja- 
cent crystals, which had registered parts of the electro- 
magnetic shower initiated by an incoming particle. But 
a crystal could only contribute to a cluster, if its energy 
deposit exceeded 2 MeV in the CB and 4 MeV in TAPS. 
The energy of the cluster was calculated by the sum of the 
energies of all contributing crystals. The cluster direction 
was determined by the weighted sum of the directions of 
the contributing crystals, using the square-root of the en- 
ergies as weight. Only clusters, which had a total energy 
of 20 MeV or higher, were used in the analysis. 

The next step in the analysis of the simulated data was 
to implement a software trigger system. This was crucial, 
since electronics and photomultipliers of the PID in the 
exit region of the CB screened TAPS from particles emerg- 
ing from the target. Some photons were converted into 
electron positron pairs in this material and subsequently 
identified as charged particles in TAPS. Protons lost parts 
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Fig. 2. The two trigger conditions. The grey shaded histograms were calculated from the simulated data. The black line 
corresponds to the experiment. Left: CB energy sum. Right: Sector multiplicity. 



of their energy in this inactive material. But these effects 
could be simulated precisely as discussed in [25] . Figure [2] 
shows the agreement between the experimental and simu- 
lated trigger conditions for selected rj — > 3tt° events. Note, 
that a sector in the trigger did not directly correspond to 
a cluster in the calorimeters. Trigger sectors had a much 
tighter restriction due to the high thresholds on single 
crystals compared to the total energy threshold on the 
clusters. Therefore, more than half of the selected events 
had a multiplicity lower than six. 

In further analysis of the experimental and the simu- 
lated data, the reaction hypothesis r\ — > 3ir° — > 67 was 
tested with a kinematic fitting technique [53]. The mea- 
sured energies E and the 9 and <j> angles were used as 
input parameters to a fit with five constraints: the invari- 
ant and the missing masses of the six identified photons 
had to be equal to the masses of the 77 meson and the 
proton, respectively, and the invariant masses of each of 
the three photon pairs had to give the tt° mass. All 15 
possible combinations to form three pairs from six pho- 
tons were tested in separate fits. Events in which at least 
one fit had a confidence level (CL) higher than 2 % were 
considered as 77 — > 3tt° events. This cut was chosen to 
reject most of the background, but lose only a small frac- 
tion of events of interest. The adjusted photon energies 
and angles from the fit with the highest CL were used to 
calculate the Dalitz plot parameter, a. 

The resolutions of the three variables E, 8 and (f> used 
in the fits were determined by a Monte Carlo simulation 
of the CB/TAPS-setup. Photons of different energies were 
generated emerging isotropically from the centre of the CB 
and the detector response was simulated in GEANT. The 
resolutions cte, erg and 0$ were then obtained by compar- 
ing the reconstructed variables with the initially generated 
values. 

Since the events were selected by cutting on the CL 
of the fits, it was important that the simulated CL distri- 
bution agreed with the measured distribution. Then cuts 
made at different confidence levels removed the same frac- 
tion of events from the measured and the simulated sam- 
ples. Figure [3] shows good agreement between the two dis- 



tributions (top left) . The rise at low confidence levels was 
produced by events where parts of the electromagnetic 
shower leaked out of the detector, such as when clusters 
were located at the edge of the CB. 

In addition, good agreement between the experiment 
and the simulation had to be achieved in the distributions 
of the variables used as constraints in the kinematic fit. 
This is also illustrated in fig [3j Here the invariant mass 
of the six photons, calculated using the directly measured 
energies and momenta, is shown for the events selected by 
the kinematic fit. There is very good agreement between 
the two distributions, which have maxima at 547.7 MeV 
and widths of a w 16 MeV. The missing mass of the mea- 
sured photon 4-momenta is shown in the lower left plot. 
Both curves have a maximum at 936.7 MeV and a width 
of a = 18.2 MeV, indicating that the chosen events satis- 
fied this constraint. The ir° mass spectrum is shown in the 
lower right plot, where the invariant masses of the three 
photon pairs from the fit with the highest CL are shown. 
The two distributions with peaks at 134.4 MeV and widths 
of a ~ 7 MeV agree within 4 %. 

As a further check of the analysis and the simulation, 
total cross sections for 77- and 7r°-photoproduction were 
determined and compared to previous experiments and 
model fits to experimental data. Figure 0] shows these 
cross sections for the energy range of 700 to 820 MeV. 
The total 77 cross section is compared to the measurement 
carried out with TAPS at MAMI [M] and the Eta-MAID 
model fit [35], which is dominated by this TAPS exper- 
iment in the threshold region. It is clearly seen that the 
results obtained agree with both, the TAPS data and the 
Eta-MAID calculation. Figure [4] (right) shows the total 
7r°-photoproduction cross section as determined using an 
analysis similar to that described above. The only dif- 
ference was that just two constraints, namely the invari- 
ant and the missing mass of two identified photons, were 
used. The 7r° cross section shows good agreement with the 
MAID2007 [3"6]l3"7] and SAID [38] models. The agreement 
is also very good in the region of the A resonance, which 
is not shown here. 
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Fig. 3. Kinematic fit results of the selected events. The grey shaded histograms and the solid lines represent the simulation 
and the experimental data, respectively. All masses were obtained from the directly measured energies and momenta (not using 
the fitted quantities). Top left: Confidence level. Top right: Invariant mass of the six photons. Bottom left: Missing mass of the 
six photons. Bottom right: Invariant masses of the three photon pairs. 
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Fig. 4. Total cross sections. Left: 7/-photoproduction cross section compared to the measurement made with TAPS at MAMI 
|34] (points with only vertical error bars) and the Eta- MAID calculation |35| (solid line). Right: 7r°-photoproduction cross section 
compared to the MAID2007 [551137] (dashed line) and SAID [38] models (solid line). 
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Although these comparisons indicate that the selected 
event sample was almost free of background, possible con- 
taminating background reactions were considered. In the 
examined energy region the main background contribution 
comes from the direct 37r°-production through the reac- 
tion 7p — > 37r°p. To study the fraction of such background 
events in the selected r/ sample 10 7 3tt° events were simu- 
lated and analysed, giving a total reconstruction efficiency 
e 37r o w 5 %. The background contribution was calculated 
using the estimate of the total cross section made in |39j , 
which resulted in <r(7p — > 37r°p) « 0.4 [ih for photon beam 
energies E 1 < 1100 MeV. It was assumed to be constant 
over the examined energy region. The contamination was 
then estimated to be 



Nrj^ 37r o a n ■ BR(r) -> 37T°) e^^ 37r o 



2%, 



(6) 



where e 3w o and e^^^o are the reconstruction efficiencies 
for the direct 37r°-production and the r\ — > 3ir° decay, 
respectively, and BR(?7 — > 37r°) is the branching ratio for 
the given decay. The total 77-photoproduction cross section 
was averaged over the observed photon beam energy range 
of 700 MeV to 820 MeV, resulting in a v » 14 fib. The 2% 
contribution is much smaller than the estimated statis- 
tical and systematic uncertainties (see section 2]) , which 
were determined to be of the order of 5 to 10% each. 
Therefore, it was neglected. Another background process, 
namely double 7r°-production with two cluster split-offs, 
was also found to be negligible. Other background pro- 
cesses were kinematically not possible due to the restricted 
tagged photon energy range of E 1 < 820 MeV. 



4 Results 
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Fig. 5. The z distribution from the r\ — > 3tt° decay events, se- 
lected by the kinematic fit. Top: Comparison of the simulated 
distributions for complete detector acceptance and 100 % de- 
tection efficiency (dashed line) and for a realistic simulation 
(solid line). Bottom: Comparison of the distributions recon- 
structed from the realistic simulation (same solid line as in the 
upper picture) and the measured data (grey shaded). 



4.1 Dalitz plot parameter 

The Dalitz plot parameter, a, for the rj — ► 37r° decay was 
determined by comparing the simulated with the mea- 
sured z distribution. The events were selected by kine- 
matic fits testing the hypothesis 7p — > ?/p — > 37r°p at 
the 2 % CL. The simulation was based on pure phase- 
space distributions with a = 0. The dashed line in fig. [5] 
(top) shows the generated z distribution for full detector 
acceptance and 100% detection efficiency. The solid line 
illustrates that realistic acceptance and efficiency already 
introduce a slope to the simulated distribution, thus, af- 
fecting our result for the Dalitz plot parameter. From the 
distribution of the difference between reconstructed and 
initially generated z values resolution for the variable z 
was found to be a z » 4.5 % and, thus, a bin width for the 
z distribution of 0.05 was chosen. 

In order to determine the Dalitz plot parameter, the 
experimental z distribution was divided by the simulated 
z distribution. This ratio R{z) was then fitted with the 
function c(l + 2az), according to eq. [3l using c and a as 
free parameters in the fit. The deviation of R{z) from pure 
phase space is then reflected in the value of a. This de- 
viation is illustrated in the lower picture of fig. [5] Here 



the simulated z distribution was normalised with c. The 
ratio R(z), also scaled with c, and the line fits for the two 
experiments are presented in fig. [6] The upper limit of the 
fit region was fixed at z = 0.9, because the last two bins 
had much poorer statistics and showed systematic devi- 
ations from a straight line. This effect occurred due to 
slightly different resolutions used in the kinematic fit for 
the experimental and the simulated data. These had to 
be applied to match the experimental and the simulated 
CL distributions. As the cut on the CL was the major re- 
striction in the presented analysis, agreement in the CL 
distributions was desired. So, the last two bins were ex- 
cluded by limiting the fit region to z < 0.9. The systematic 
effect of these two bins on the Dalitz plot parameter was 
examined in a test, fitting a line to the region < z < 1 
(see section |4~2"|) . 

The two experiments differed in the CB energy sum 
trigger threshold of « 390 MeV and £ t 2 hr w 60 MeV, 
and the tagged photon energies, 680MeV< <820MeV 
compared with 200MeV< < 820 MeV. The results 

of the standard analysis described above for these two 
experiments, ct\ and 0:2, are listed in the first line of table 
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Fig. 6. Ratio R(z) of the simulated and the measured z distri- 
butions for the two experiments. The Dalitz plot parameter, a, 
was obtained directly from the fit of the function c(l + 2az). 
In the figure, R(z) has been scaled with c so that the fitted 
lines have intercepts equal to 1 on the _R(z)-axis. Top: For the 
rj decays experiment. Bottom: For the radiative 7r°-production 
experiment. 



[H The experimental statistics and the reduced x 2 -values 
(x 2 /ndf) derived from the fits to the z distribution are also 
given. The listed statistical uncertainties were taken from 
the errors of the fits. The results from the two experiments 
agree very well. 



4.2 Systematic uncertainty 

The systematic uncertainties were estimated in a series 
of tests, varying only one analysis parameter at a time. 
The results from all tests for both experiments are listed 
in table El For each test the difference compared to the 
standard analysis is indicated. 

The first test demanded detection of exactly one pro- 
ton in addition to the six photons as an event selection cri- 
terion (67+lp), requiring accepted events to have a detec- 
tor cluster structure consistent with the 7p — * rjp — * 37r°p 
reaction. The results of this test show, within one stan- 
dard deviation, good agreement with the standard values 
for both experiments. 

Then there was a series of tests where just the cut 
on the CL was altered between 1% and 50%. For the 



standard analysis it was set to 2 %. With the cut on the 
lower value of I % the effect of relaxing the acceptance 
criteria, and possibly including some events from back- 
ground processes, was investigated. It showed small devi- 
ations compared to the standard analysis. Raising the cut 
to 5%, 10%, 20% and 50% led to a higher purity of the 
chosen event sample. Most results show slight differences 
compared to the standard value. Cuts on high confidence 
levels led to larger deviation. But, since at the same time 
the statistical significances decreased, the results are still 
within one standard deviation compared to the standard 
analysis. 

The trigger threshold for the CB energy sum was de- 
termined for one experiment to be E^ hl « 390 MeV. Since 
this value could not be fixed exactly, in a series of tests it 
was varied around the standard value. But these tests were 
only performed for the experiment with the higher cut, 
because the threshold of the second experiment, roughly 
60 MeV, was so low that almost no rj — > 37r° events were 
rejected (compare fig. [2]). The results of these tests showed 
good agreement with the standard value. 

In another test TAPS clusters were ignored in the anal- 
ysis (no TAPS), so that all six required photons had to be 
detected by the CB. This was a check, if the acceptance in 
forward directions was simulated precisely enough. Such 
a test was necessary, because a lot of inactive material in 
the exit region of the CB screened TAPS from particles 
emerging from the target. Both results are slightly lower 
than the standard values, but are within one a compared 
to them. 

Possible systematic effects from misidentification of 
photons as charged particles were studied in an analy- 
sis (no ID) omitting the particle identification methods 
described above. All clusters, registered in the CB and 
TAPS, were accepted as photon candidates. Practically 
no deviations from the standard values were found. 

As mentioned above, in a test the fit region was ex- 
tended to z = 1.0. With the fit over the full z range the 
influence of the last two bins at high z values was investi- 
gated. The results in table[3]show no effect for the second, 
but a drop in a for the first experiment. Also the reduced 
X 2 indicates the the description of a straight line does not 
hold here for the full z range. But within one standard de- 
viation the test results are compatible with the standard 
values. 

The upper limit of z = 0.75 was chosen for a test to 
investigate possible influences of the region in the Dalitz 
plot with z > 0.756, where statistics significantly decrease 
(see fig. [SJ . In eq. 0] z is given in such a way that all points 
with z < 0.756 have the same probability for the pure 
phase space process. In [40] it is shown that at z = 0.756 
a cusp in R{z) arises. The result of this test indicates just 
minor changes compared to the standard values. 

As discussed in [40], the ratio R{z) exhibits two more 
cusps, one at z = 0.597, the minimum value to reach 
the 7r + 7r~ threshold in the Dalitz plot, and another at 
z — 0.882, corresponding to the maximum value to touch 
this line. These cusps arise due to the ir + ir~ — > 7r°7r° 
charge exchange-reaction, which also produces a cusp in 
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Table 3. Results from tests concerning the systematic uncertainty of the extraction of the Dalitz plot parameter, a, for the 
rj — > 37T° decay. Given are abbreviations for the test conditions, which are explained in the main text, the test results with their 
statistics N and the \ 2 /ndf of the straight line fit for both experiments. 



Test 




m /in 3 

ivl / 1U 


X / ndt 


Ot-2 


\r„ /in 3 

iV2 / 1U 


X /nat 


1 


Standard 


-0.0315 ± 0.0020 


1120 


20/16 


-0.0324 ± 0.0024 


724 


10/16 


2 


u lily 


-0.0323 ± 0.0033 


430 


15/16 


-0.0357 ± 0.0035 


322 


18/16 


3 


CL=1 


-0.0308 ± 0.0020 


1177 


20/16 


-0.0331 ±0.0023 


759 


9/16 


4 


CL=5 


—0.0327 ± 0.0021 


1022 


22/16 


-0.0336 ± 0.0025 


663 


12/16 


5 


CL=10 


-0.0328 ± 0.0022 


919 


21/16 


-0.0322 ± 0.0026 


600 


11/16 


6 


CL=20 


-0.0332 ± 0.0024 


774 


17/16 


-0.0296 ± 0.0028 


510 


16/16 


7 


CL=50 


-0.0295 ± 0.0031 


453 


20/16 


-0.0297 ± 0.0036 


303 


11/16 


8 


E thl = 380 MeV 


-0.0323 ± 0.0020 


1121 


20/16 








9 


E thl = 388 MeV 


-0.0316 ± 0.0020 


1120 


20/16 








10 


E thl = 392 MeV 


-0.0313 ± 0.0020 


1119 


20/16 








11 


E thl = 400 MeV 


-0.0307 ± 0.0020 


1117 


21/16 








12 


no TAPS 


-0.0300 ± 0.0024 


797 


13/16 


-0.0316 ± 0.0028 


516 


11/16 


13 


no ID 


-0.0311 ± 0.0023 


838 


12/16 


-0.0306 ±0.0027 


556 


12/16 


14 


2 < 1 


-0.0295 ± 0.0019 


1149 


31/18 


-0.0322 ± 0.0022 


742 


13/18 


15 


z < 0.75 


-0.0319 ± 0.0024 


1005 


19/13 


-0.0312 ±0.0029 


650 


8/13 


16 


z < 0.6 


-0.0284 ± 0.0034 


812 


13/10 


-0.0277 ±0.0041 


524 


7/10 



the 7r°7r° mass distribution of the r\ — > 3tt° decay 001,42, 
|4"3"1|44U45U4"6] . The cusp at higher z in the ratio R{z) was 
already excluded by the standard analysis. Reduction to 
z = 0.6 makes a substantial decrease in the result for a, 
but the statistical error is also larger. In [3] Ditsche, Kubis 
and MeiBner calculated that, if the Dalitz plot parameter 
is obtained by a fit over the range < z < 0.597, a drop of 
5 % in the absolute value of a should be visible, compared 
to a fit over the full z range. The results from test 16 show 
a decrease of roughly 10% and 15% for the two experi- 
ments. But within the increased errors, the test results 
still agree with the standard values. So, no reliable state- 
ment on the influence of the cusps in R(z) can be made 
and the differences in a were just taken into account in 
the determination of the systematic uncertainties. 

The systematic uncertainties were then calculated for 
both experiments separately according to 



^syst(ai) 



Efc Kfc ~ a i) ■ N H 



(7) 



where on stands for the standard values of the two exper- 
iments. The a' ik are the results of the different tests listed 
in table [3] with their statistics Nik ■ Thus, the systematic 
uncertainties were calculated by the sum of the differences 
between the test results and the standard values, weighted 
with the statistics of the tests. Positive and negative de- 
viations were handled separately in this process. 



4.3 Final results 

As final results 

OL2 ■ 



-0.0315 ±0.0020 



-0.0012 
-0.0009 



-0.0324 ±0.00241°™^ 



(8) 
(9) 



were obtained. The result of the analysis presented in this 
paper was then calculated by the weighted mean of the two 
values, while the inverse variances were taken as weights. 
The statistical uncertainty was calculated from 



1 



E,(iM 2 



(10) 



The systematic uncertainty was taken from the highest 
systematic uncertainty given above. As final result we 
quote 

a = -0.032 ± 0.002 stat ± 0.002 syst (11) 

This value is consistent with the current PDG average 
[10j . which is dominated by the Crystal Ball result mea- 
sured at BNL, and agrees reasonably well with the results 
published by the KLOE [15], CELSIUS/WASA [20] and 
WASA-at-COSY [21] collaborations. 

After the upgrade of MAMI to 1.5 GeV maximum elec- 
tron energy [47] , a new experiment on the neutral decays 
of the rj meson was performed in the tagged photon energy 
range of 700 MeV to 1400 MeV. An independent analysis 
of this new experiment [48 obtained approximately 3 • 10 6 
rj —* 37T° events, giving a value for a, which is in very good 
agreement with the result presented in this paper. 



5 Summary 

The Dalitz plot parameter, a, for the rj — ► Sir decay has 
been measured with the Crystal Ball detector and TAPS 
at the electron accelerator facility MAMI-B in Mainz. The 
•q mesons were produced by bremsstrahlung photons emit- 
ted by the 883 MeV electrons and detected by the Glasgow 
tagged photon spectrometer. The kinematic fitting analy- 
sis selected 1.8- 10 6 7p — > rjp — > 37r°p events and the result 
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a — —0.032 ± 0.002 sta t ± 0.002 syst is in good agreement 
with other high- precision experiments. A possible influ- 
ence on the Dalitz plot parameter by the 7r + 7r~ — > 7r°7r° 
contribution to the amplitude was found to be small, but 
is included in the systematic uncertainty. 
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